Spin Torque in Anisotropic Tunneling Junctions 
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Spin transport in magnetic tunnel junctions comprising a single magnetic layer in the presence of 
interfacial spin-orbit interaction (SOI) is investigated theoretically. Due to the presence of interfacial 
SOI, a current-driven spin torque can be generated at the second order in SOI, even in the absence 
of an external spin polarizer. This torque possesses two components, in-plane and perpendicular to 
the plane of rotation, that can induce either current-driven magnetization switching from in-plane to 
out-of-plane configuration or magnetization precessions, similarly to Spin Transfer Torque in spin- 
valves. Consequently, it appears that it is possible to control the magnetization steady state and 
dynamics by either varying the bias voltage or electrically modifying the SOI at the interface. 



PACS numbers: 72.25.-b,73.43. Jn,73.40.Rw,73.43.Qt 

Since the theoretical prediction and experimental 
observation of current-driven magnetization control^, 
spintronics is considered as a promising candidate for 
low energy consumption devices^. Up until now, Spin 
Transfer Torque (STT) has essentially been observed 
in inhomogeneous magnetic structures such as Spin- 
Valves, Magnetic Tunnel Junctions and Magnetic Do- 
main WallsA. Moreover, recent experimental obser- 
vations on single ferromagnet-bascd structures have 
opened promising opportunities for new device concepts. 
For example, Tunneling Anisotropic Magnctorcsistancc 
(TAMR) has been observed in Magnetic Tunnel Junc- 
tions (MTJs) comprising a single ferromagnet^— (re- 
ferred to as Semi-MTJs - SMTJs). Ref. i proposed to 
exploit the spin-orbit interaction (SOI) present in a single 
ferromagnetic layer to electrically control the magnetiza- 
tion directio n 10 ' 11 . Alternatively, the voltage-controlled 
manipulation of magnetic anisotropy of thin magnetic 
layers through thick insulators has been achieved^. Effi- 
cient electrical control of the magnetization direction of a 
single ferromagnetic layer combined with sizable TAMR 
effect would offer powerful perspectives for spin-based 
memory devices. 

The key ingredient of TAMR£~£, SOI-induced torque 
(SOI-ST)2 and magnetic anisotropyi^ is the SOI aris- 
ing either from Bulk Inversion Asymmetry^! (BIA) or 
Structure Inversion Asymmetr y 10 ' 12 (SIA). In the latter 
case, the presence of a large potential gradient at the 
interface between, say, an insulator and a metal, gener- 
ates a local electric field perpendicular to the interface 
W = — Ez, inducing a SOI of the form^ 



Hft(k) = Bfl(k) • <r = aff(k X z) ■ &. 



(1) 



The precise form and magnitude of an terms has been 
widely studied in semiconductor 2DE G 13 ' 14 , and only re- 
cently at metallic interfaces and can be as large as 1- 
5eV.A 2 ^. 

In this Letter, we suggest that the presence of such 
an interfacial SOI at the interface between a ferromagnct 
and an insulator in a SMTJ is responsible for a non- 
equilibrium spin torque. We demonstrate that the mag- 
netization of the ferromagnetic layer of such a junction 



can be controlled or excited by an external bias volt- 
age applied across the junction. In the non-equilibrium 
regime, the interfacial SOI generates a bias voltage- 
driven spin torque on the ferromagnetic layer inducing 
either magnetization switching or self-sustained magnetic 
precessions. This effect belongs to the family of SOI- 
induced spin torques*^ and the correspondences with 
the conventional STT observed in spin-valves^ will be 
discussed at the end of this paper. 
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FIG. 1. (Color online) Left panel: [top] Schematics of a Semi- 
Magnetic Tunnel Junction, [bottom] Influence of the spin 
torque on the magnetization; Right panel: Angular depen- 
dence of (a) TAMR and (b) perpendicular torque T± at zero 
bias, for a rotation in the (010) planes and olr £[l,5]eV.A 2 . 
The parameters are adapted for Fe/MgO interface: C/o=leV, 
d=0.6nm, fct = 1.09nm -1 , k l F = 0.4nnT\ Inset: <£- 
dependence of the perpendicular torque. 



In regular MTJs, the transport is essentially con- 
trolled by the interfacial densities of states^!. There- 
fore, the tunneling magnetoresistance (TMR) and STT 
are proportional to the interfacial polarization of the 
electrodes, P,: TMR cx PlPr, Tii oc P^sin^, and 
ti oc PlPr sin flii. In contrast, in the case of SMTJs, 
P L = and both STT and TMR should vanish. How- 
ever, in the presence of interfacial SOI it has been found 
that the resistance of the junction depends on the orien- 
tation of the magnetization against the interface^. 
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This TAMR effect has been observed using 
semiconductors^, metals^ and hybrid structures^. The- 
oretical investigations^^ have shown that TAMR is well 
described with Rashba and Dresselhaus effective SO I 13 i 19 
(the latter being restricted to non-centrosymmetric crys- 
tals). Following Eq. (p}, Rashba SOI induces an angular 
momentum transfer from the spin momentum to the 
orbital degree of freedom. Since the interfacial SOI is 
linear in k, this effect vanishes at the first order in SOI 
after averaging over the Fermi sphere [(H R (k)) = 0]. 
But in a ferromagnet, due to the presence of the local 
magnetic exchange, a net transfer of angular momentum 
occurs at the second order [(iJij(k) 2 ) ^ 0]. As a 
consequence, magnetic anisotropy as well as TAMR 
appear, proportional to the second order in Rashba SOI, 
oc a ^ 18 ! 19 . Interestingly, since the angular momentum 
transfer involves only in-plane components of the spin 
and orbital momentum (i.e. o~ x . y and k x ^ y ), the itinerant 
spin density produced by this transfer lies in the (x,y) 
plane only. 

In the present work, we consider the F/I/N trilayer 
depicted in Fig. [TJ where F is a ferromagnetic layer 
(Co, Fe, Ni and compounds), I is an insulator (MgO, 
AlOx, GaAs) and N is a normal metal (Cu, Ag, Au, Pt 
etc.) . In order to capture the most relevant features of the 
mechanism described here, we choose a minimal model 
only considering the most pertinent material parameters. 
Matos-Abiague ct alJ£ showed that in the case of cen- 
trosymmetric barriers (such as AlOx, MgO), the TAMR 
is mostly due to the Rashba SOI at the interface between 
the ferromagnet and the tunnel barrier. The free electron 
Hamiltonian of the junction is then 

h 2 1 

H = + U{z) + H R 5{z- z L ). (2) 
2 m(z) 

Hr is given in Eq. ((T|), m(z) is the effective mass of the 
electron, equal to m in the electrodes and m e ffmQ in 
the barrier. U(z) is the potential of the junction, given 
by 

U, <0 = J<x-M + ^, U z>d = -^ (3) 

U 0<z<d = U + (i - ^)eV b , (4) 

where Uq and d are the barrier height and thickness, 
Vb is the bias voltage, J is the s — d exchange cou- 
pling, <x is the vector of Pauli spin matrices and M = 
(sin 9 cos 4>, sin 9 sin <j>, cos 9) is the magnetization direc- 
tion of the ferromagnetic electrode [see Fig. [Q. 

The mechanism giving rise to itinerant spin density can 
be understood by looking at the spin density continuity 
equation ^ = -h([&,H]), which reads 

dm ih_ 2 J 1 ~ ... 

— = — V • <r ® V) - — m x M + — ( [tr,H R .(5) 
dt m n in 

The first term is the regular spin current divergence in the 
absence of SOI, the second term is the torque exerted by 



the itinerant spin density m on the local magnetization 
M and the last term is the torque between the itinerant 
spin and orbital angular momentum. For instance, in 
MTJs the last term is generally zero and the spin torque 
is directly associated with the spatial variation of the 
spin current cx (<r <£> V^£. In the present case, Rashba 
SOI acts like a source for spin density and therefore the 
spin torque is no more simply related to the spin current. 
Consequently, the proper way to evaluate the spin torque 
is to calculate directly the local itinerant spin density 
m(z). We assume a semi-infinite magnetic layer, as is 
usually done in magnetic tunnel junctions^: The total 
spin torque exerted by the transverse spin density m on 
the magnetization M is defined as T = J J v m x 'M.dV, 
where V is the volume of the magnetic layer. 

The charge and spin currents are then evaluated using 
the conventional definitions 

Je = i 3 E / ^kii^r^kilR^^kn)/,:], (6) 

s,i 

m = E / dEd^riE^^m^Ofi, (7) 

*i>f (E, k| | ) being the Hartree-Fock two-component wave 
function for an electron of energy E and in-plane wave 
vector k||, issued from the i-th reservoir with a Fermi- 
Dirac distribution /j. Assuming convenient material 
parameters (barrier characteristics and effective mass), 
the present minimal model satisfyingly reproduces the 
TAMR obtained by Matos-Abiague et al&i^ in the case 
of Fe/GaAs/Au and the one obtained by Park et alii 
in (Pt/Co) n /A10x/Pt, with a reasonable degree of ac- 
curacy. The angular dependence of the zero bias con- 
ductance for rotation of the magnetization in the (010) 
plane for an = 1 — 5eV.A 2 is displayed in Fig. HJa) and 
shows the expected (cos 29 — 1 )-dcpendcnc o 6 i 18 . This an- 
gular dependence can be understood by noticing that the 
junction is physically equivalent upon the transformation 
6 -> -9 9 + 7ri£. 

Since the torque is by definition perpendicular to the 
magnetization, it reads 

T = Tj|M x (z x M) + T±z x M. (8) 

The spin torque possesses two components that can be 
referred to as in-plane 21 1 , and perpendicular torques, T±. 
At zero bias, only the perpendicular torque T± is non- 
zero and displays an angular dependence on the form 
sin 29 [Fig. H[b)] which favors the perpendicular config- 
uration. As mentioned above, the itinerant spin density 
due to Rashba SOI lies in the (x,y) plane of the junction 
[Eq. (pQ)]. Consequently (i) when the magnetization is 
oriented along z (9 = 0), no transfer occur ((o- x ,y) = 0) 
and the spin torque vanishes; (ii) when the magnetiza- 
tion is oriented in the (x,y) plane (9 = tt/2), since the 
spin density lies itself in this plane, the spin torque is also 
zero. Finally, the junction is invariant under 0-rotation, 
so the spin torque (and resistance) does not depend on 
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[Fig. QJb), inset]. This gives rise to an angular de- 
pendence on the form sin 29 as well as four associated 
stable magnetic states: two perpendicular to the plane 
(9 = 0, 7r) and two in the plane (9 = ±7r/2). 
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FIG. 2. (Color online) Non-equilibrium in-plane torque 
Tj|(Vt) as a function of (a) Rashba parameter and (b) bias 
voltage. The dashed line shows the non-equilibrium perpen- 
dicular torque for V b =0.8V (a) and a R = 4eV.A 2 (b). The 
parameters are the same as in Fig. [TJ 

When applying a bias voltage across the junction, spin 
polarized electrons tunnel through the barrier, signifi- 
cantly modifying the spin imbalance in the left electrode. 
Therefore, both the TAMR and the magnetic anisotropy 
are strongly affected by the voltage. In Fig. [2] the 
non-equilibrium in-plane torque, Tj|(Vf,), is represented 
a function of an [Fig. |2fa)] for different bias voltages, 
and as a function of the voltage [Fig. Hfb)] for different 
aR. The spin torque displays a quadratic dependence on 
the Rashba parameter as expected from the symmetry of 
the Rashba Hamiltonian (linear in k). For larger Rashba 
parameters, higher orders in an appear (not shown). In- 
terestingly, the amplitude of the non-equilibrium perpen- 
dicular torque T±(Vb) — T±(0) is about one order of mag- 
nitude smaller than the in-plane torque [dashed line in 
Fig. HJa) and (b)]. 

The form of the torque displayed in Eq. © is similar to 
the usual STT in a MTJ whose polarizer is oriented along 
ail. Whereas the perpendicular torque oc z x M com- 
petes with the demagnetizing field and the perpendicular 
anisotropy, the in-plane torque «Mx(zxM) competes 
with the damping. As a consequence, one expects the 
current-driven torque to produce current-induced mag- 
netization switching from out-of-plane to in-plane and 
vice-versa as well as current-driven magnetization pre- 
cessions. 

In order to illustrate the current-driven magnetiza- 
tion dynamics that can be generated in such structures, 
we numerically solve the macrospin Landau-Lifshitz- 
Gilbcrt equation in the presence of spin torque using 
a fourth-order Runge-Kutta method^. A ferromagnetic 
layer is considered, possessing both in-planc and out-of- 



> 

a. 

o 
i_ 
+j 
o 

■i2 0.15 
c 

< 

0) 

c 

ro 

cl 
c 



0.10 



o 

N 

E 



0.05 



0.00 



Precessions 





inn \OPP 

IPP valM 




V I - 3 


Stable 


\ \ 
\ \ 


Out-ot-Plane 


\ \ 




\ \ 




\ \ 




Stable N V 




In-Plane \* 





O -0.03 -0.02 -0.01 0.00 0.01 0.02 0.03 
Normalized In-Plane Torque 

FIG. 3. (Color online) Stability Diagram of the magnetization 
in the presence of perpendicular and in-plane anisotropy, as 
well as in-plane torque. The parameters are magnetic damp- 
ing rj = 0.01, perpendicular anisotropy Q± — 0.75, saturation 
magnetization M s — 1.6T— . The dashed lines are guides for 
the eye. 



plane magnetic anisotropies Q| = Hff JM S and Q± = 
/M s , where Hff (H^) is the in-plane (perpendicular) 
anisotropy field and M s is the saturation magnetization. 
Both anisotropies are needed in order to achieve stable 
in-plane and out-of-plane magnetization states. Since 
Tj_ < < T||, the non-equilibrium perpendicular torque will 
be disregarded. In addition, no external field is applied 
and only the in-plane torque Tj| = Tji/ (figM^d) is con- 
sidered. Finally, for the numerical simulations, we adopt 
parameters close to the one measured by Ikeda et ali2£ 
(see Fig. EJ). 

Fig. [3] displays the normalized resistance of the junc- 
tion (oc cos 29) taken at a time t=35ns for a magneti- 
zation initially perpendicular to the plane. Similar re- 
sults are found when the magnetization is initially in- 
plane (not shown). Interestingly, four zones can be dis- 
tinguished. For large positive in-plane torque, Tj| > 
T(h = 77(1 + Q|| — Q±) [rj is the Gilbert damping), the 
magnetization is stable in the perpendicular direction. 
When T|| < Tth, out-of-plane [OPP - Fig. 0Ja)] and in- 
plane [IPP - Fig. H[b)] magnetization precessions appear, 
while for even larger negative in-planc torque the in-planc 
anisotropy overcomes the perpendicular anisotropy and 
in-plane stable states are reached [see Fig. Hfc)]. The 
dependence of the resonance frequency as a function of 
the in- plane torque is reported on Fig. EJd), where the 
positive slope (left hand) is attributed to out-of-plane 
precessions, whereas the negative slope (right hand) is 
attributed to in-plane precessions. 

The present current-driven spin torque can be read- 
ily compared with the conventional STT— and the SOI- 
ST— First, the present torque has its origin in the 
interfacial SOI, rather than in the inhomogeneous or dis- 
continuous magnetic texture so that no external polar- 
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FIG. 4. (Color online) Magnetization trajectory for a mag- 
netization initially perpendicular to the plane and submitted 
to an in-plane anisotropy Qn = 0.2 and (a) tii = —0.0125, 
(b) r| = -0.02, (c) t n = -0.0375; (d) Dependence of the 
oscillation frequency as a function of the in-plane torque. 



produces only an effective magnetic field. Furthermore, 
similarly to the conventional STT, the in-plane torque T| 
competes with the damping and can excite self-sustained 
magnetic precessions in the absence of external magnetic 
field, provided that stable in-plane and out-of-plane mag- 
netization stable state can be achieved. 

For the parameters exploited here, the spin torque can 
be as large as 100eV/^ m 2 (sec Fig. 2(b)), which is 
comparable to spin transfer torque in MTJ$i£. From 
an experimental point of view, the key element needed 
to observe such spin torque is a large spin-orbit split- 
ting at the interface between the ferromagnet and the 
barrier. As mentioned in the introduction, a number of 
results have been produced on SMTJs based on metals^ 
or semiconductors^. A good rule of thumb to design such 
efficient interfaces is to maximize the TAMR, as shown in 
Ref.— . An efficient procedure to detect such a torque in 
the case of small bias (or small Rashba parameter) is to 
investigate spin torque-driven ferromagnetic resonance 2 ^ 
or spin-diode effects^. 
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